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ABSTRACT 

We observed the first known very high energy (VHE) 7-ray emitting unidentified source, TeV 
J2032+4130, for 94 hours with the MAGIC telescope. The source was detected with a significance of 
5.6cr. The flux, position, and angular extension are compatible with the previous ones measured by the 
HEGRA telescope system five years ago. The integral flux amounts to (4.5±0.3stat±0.35sys) x 10~^^ 
ph cm~^ s~^ above 1 TeV. The source energy spectrum, obtained with the lowest energy threshold 
to date, is compatible with a single power law with a hard photon index of r=-2.0±0.3stati0.2sys. 
Subject headings: gamma rays: observations 



1. INTRODUCTION 

The TeV source J2032+4130 (Aharonian et al. 2002) 
was the first unidentified very high energy (VHE) 7-ray 
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source, and also the first discovered extended TeV source, 
hkely to be Galactic. 

Intensive observational campaigns at different wave- 
lengths have been carried out on TeV J2032-I-4130. Butt 
et al. (2003) presented an analysis of the CO, HI, and 
infrared emissions, together with first observations by 
Chandra (5 ksec) and a reanalysis of VLA data. These 
observations showed that the TeV source region is posi- 
tionally coincident with an outlying group of stars (from 
the Cygnus 0B2 core), although they failed to identify a 
counterpart. Mukherjee et al. (2003) analyzed the same 
Chandra data and provided optical follow-up observa- 
tions of several of the brightest X-ray sources, confirm- 
ing that most were either O stars or foreground late-type 
stars. A deeper Chandra observation (50 ksec. Butt et 
al. 2006), found hundreds of star-like sources and yet no 
diffuse X-ray counterpart emission. 

A deep 50 ksec) XMM-Newton exposure has also 
been obtained (Horns et al. 2007). After the subtrac- 
tion of the contribution of known sources from the data, 
an extended X-ray emission region with a FWHM size 
of ~ 12 arcmin was reported. The centroid of the emis- 
sion is co-located with the position of TeV J2032-I-4130 
and was proposed as the counterpart of the TeV source. 
The question whether the result reported by Horns et al. 
can be interpreted as a truly diffuse background, or it 
could be a result of unresolved X-ray sources, remains 
disputable. 

Paredes et al. (2007) and Marti et al. (2007) 
have provided deep radio observations covering the TeV 
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J2032+4130 vicinity using the Giant Metrewave Radio 
Telescope and discovered a population of radio sources, 
some in coincidence with X-ray detections by Butt et al. 
(2006) and with optical/IR counterparts. At least three 
of these sources are non-thermal, and one has a hard X- 
ray energy spectrum. They found extended non-thermal 
diffuse emission in the radio band apparently connecting 
with one or two radio sources. It is yet to be determined 
if one or more of these sources is similar to some of the 
known 7-ray binaries (e.g., Aharonian et al. 2006, Albert 
et al. 2006a). 

Several theoretical explanations for the TeV emission 
from J2032-I-4130 have been given. Among them, those 
related with extragalactic counterparts, e.g., a radio- 
galaxy (Butt et al. 2006) or a proton blazar (Mukherjee 
et al. 2003) , face the difficulty of explaining the extended 
appearance of the source. Gamma-ray production in hy- 
pothetical jet termination lobes of Cyg X-3 was explored 
(Aharonian et al. 2002), but the putative northern lobe 
of Cyg X-3 (now considered a mere thermal HII region, 
Martf et al. 2006) is far from the location of the TeV 
source. A yet unknown pulsar wind nebula (PWN) was 
proposed by Bednarek (2003), although no clear PWN 
signal was observed. A distant microquasar was pro- 
posed by Paredes et al. (2007), perhaps related with 
one of the X-ray/radio sources they discovered. If such 
an association is accepted, the extension of the source 
could be explained by the diffusion of accelerated parti- 
cles into a hypothetical nearby molecular enhancement 
(see Bosch- Ramon et al. 2005). Torres et al. (2004) 
and Domingo-Santamarfa & Torres (2006) studied the 
relationship between the TeV emission and the known 
massive stars in the area, through the interaction of rel- 
ativistic protons with wind ions. The distribution of stars 
in the neighborhood favors this interpretation (Butt et 
al. 2006). An explanation involving the excitation of 
giant dipole resonances of relativistic heavy nuclei in ra- 
diation dominated environments has also been suggested 
(Anchordoqui et al. 2007). 

2. PREVIOUS VERY HIGH ENERGY 7-RAY 
OBSERVATIONS 

We start by making a brief summary of what has been 
claimed by other experiments observing at the highest 
energies. 

The HEGRA lACT using four years of data (from 
1999 to 2002) found a source to the north of Cygnus 
X-3, steady in flux over the years, extended, with radius 
6.2±1.2stat±0.9sys arcmin, and exhibiting a hard energy 
spectrum with a photon index of F— -1.9± 0.1stati0.3sys 
(Aharonian et al. 2005). Its integral flux above 1 TeV 
amounts to ~5% of the Crab Nebula, assuming a Gaus- 
sian profile for the intrinsic source morphology. The cen- 
ter of the source position was determined quite accu- 
rately at aj2ooo=20'^ 31™ 5T.0 ± 6!2stat ± 13';7sys and 

'5j2000=41°29'56'.'8± l.'lstat ± l%ys- 

The Whipple collaboration reported an excess at the 
position of the HEGRA unidentified source (3.3f7) in 
their archival data of 1989 and 1990 (Lang et al. 2004), 
with a flux level of ~ 12% of the Crab Nebula for E>600 
GeV. The detected flux is in conflict with the HEGRA 
flux level and steady nature of the source, assuming 
they all have the same origin. This large difference be- 
tween the detected flux levels, if physical, might suggest 



episodic emission (with low duty cycle) or variability over 
timescale measured in years. Nevertheless, the existence 
of 7-ray variability is difficult to reconcile with the ex- 
tended appearance of the source. Also the large differ- 
ence might be in part due to unspecified systematic errors 
on the fiux determination. Recently, the Whipple collab- 
oration reported new observations of this field done with 
their 10-m telescope for 65.5 hours during 2003 and 2005 
(Konopelko et al. 2007). Their data is consistent with 
either a point-like or an extended source with less than 
6' angular size. Regarding the position, the HEGRA and 
the latest Whipple data are barely in agreement: Their 
centers of gravity are ^9' apart, and only agree when 
adding up the spatial uncertainties in both data sets in 
opposite directions. Konopelko et al. do not provide a 
energy spectrum for this source, but give a 8% Crab-level 
fiux (although with no energy threshold specified) under 
the assumption of a steep (Crab-like) energy spectrum. 

The Cygnus region shows an excess in the Milagro 
data (Abdo et al. 2006). The flux at 20 TeV in a 3x3 
square degree region centered at the HEGRA position is 
(9.8±2.9stat±2.7sys) X 10^15 TeV^^ cm^^ g-i assuming 
a differential energy spectrum E^^-^. This flux is three 
times the HEGRA fiux extrapolated at 20 TeV. The Ti- 
bet Air Shower detector recently reported evidence for an 
excess also in their VHE 7-ray candidate set from this 
region (Amenomori et al. 2006). 

In this rich observational and theoretical context we 
report here on MAGIC telescope observations of TeV 
J2032-I-4130. 

3. MAGIC OBSERVATIONS AND RESULTS 

The MAGIC single dish Imaging Air Cherenkov Tele- 
scope (see e.g., Cortina et al. 2005 for a detailed descrip- 
tion) is located on the Canary Island of La Palma. Its an- 
gular (energy) resolution is approximately 0.09° (20%), 
and the trigger (analysis) threshold is 55 (60) GeV at 
zenith in dark conditions (see Albert et al. 2007a). One 
of the unique characteristics of MAGIC is its capability of 
observing under moderate Moon light illumination (Al- 
bert et al. 2007b) albeit with a slightly elevated thresh- 
old. 

The field of view of TeV J2032+4130 was observed 
with MAGIC for more than 100 hours distributed in 
2005, 2006 and 2007, see table □ During the first pe- 
riod in Summer 2005, the observation was carried out 
in ON/OFF mode, that is, the source was observed on- 
axis while observations from an empty, nearby field of 
view were used to estimate the background. In Summer 
2006 and 2007, the data were taken in Wobble mode, us- 
ing five positions around the HEGRA position instead of 
the usual two symmetrical position in order to monitor a 
wider field of view. Quality cuts based on the trigger and 
after-cleaning rates were applied in order to remove bad 
weather runs and data spoiled by car or satellite light 
flashes. After these quality cuts the total observation 
time is 93.7 h. The energy range for which we report 
these results is signiflcantly above the aforementioned 
trigger and analysis threshold energies due to the fact 
that the observations were scheduled during moonlight 
and at relatively high zenith angles (up to 44°). 

The data analysis was carried out using the standard 
MAGIC analysis and reconstruction software (Bretz & 
Wagner 2003). It follows the general stream explained 
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in Albert et al. (2006b, c,d). After calibration and two 
levels of image cleaning tail cuts (for image core and 
boundary pixels, see Fegan 1997), the camera images are 
parameterized by the so-called image parameters (Hillas 
1985). The Random Forest method was applied for the 
7/hadron separation (Albert et al. 2007c). Using this 
method a parameter, dubbed hadronness (H), can be 
calculated for every event and which is a measure of 
the probability that the event is not 7 like. The 7 like 
sample is selected for images with a H below a specified 
value, which is optimized using a sample of Crab Neb- 
ula data processed with the same analysis stream. An 
independent sample of Monte Carlo 7-showers was used 
to determine the cut efficiency. Since part of our obser- 
vations was recorded during partial moon-shine, we have 
corrected the efficiency loss due to the increase of ambi- 
ent light following the procedure outlined in Albert et al. 
(2007b). 

The ^^-distribution was calculated, being 6 the angu- 
lar distance between the source direction and the recon- 
structed arrival direction of the showers. The reconstruc- 
tion of individual 7-ray arrival directions makes use of 
the DISP method (Domingo-Santamaria et al. 2005). 
The expected number of background events is calculated 
using five regions symmetrically placed for each wobble 
position with respect to the center of the camera and ref- 
ered to as anti-sources. Figure[l]shows the distribution of 
the 9^ parameter for the excess observed from the direc- 
tion of the source, for a SIZE cut of 800 photoelectrons 
(pe) . This relatively high SIZE cut was selected in order 
to optimize the sensitivity for a source with such a hard 
energy spectrum observed during moonlight. Therefore, 
the total number of 7-like excesses after Hillas cuts and 
applying a cut in 9'^ <0.05, is A'cx = 233, for which a 
total significance of 5.6(t is obtained. Table [5] shows the 
number of excesses above background for the different 
observing periods. 

The excess is fitted to a Gaussian function folded with 
the telescope PSF, as obtained from Monte Carlo sim- 
ulations and validated with Crab Nebula observations. 
The source is extended with respect to the MAGIC 
PSF. Its intrinsic size assuming a Gaussian profile is 
Csrc =5.0±1.7sta±0.6sys arcmin. The exact shape of 
the source, even if similar to the keV diffuse emission 
reported by Horns et al. 2007, cannot be completely 
trusted due to limited statistics and telescope pointing 
systematics. 

Figure [2] shows the Gaussian-smoothed ((7=4') map 
(0.65° X 0.65°) of 7-ray (background subtracted) around 
TeV J2032-F4130 for energies E > 500 GeV. The po- 
sition of a few previously observed 7-ray source can- 
didates are also shown, namely Cyg X-3, the EGRET 
source 3EG J2033-I-4118 (with its confidence contour at 
95%), the Wolf Rayet star WR 146, and the Whipple 
and HEGRA experimental positions. The regions around 
Cyg X-3, WR 146 and 3EG J2033+4118 have been fur- 
ther investigated by us and no detection is obtained for 
a steady emission. The upper limit fluxes (Rolke et al. 
2005) for 95% confidence level, above 500 GeV for a 
point-like source at these positions are given in table [3] 

To determine the best position of the MAGIC detec- 
tion the excess map was fitted to a 2D bell-shaped func- 
tion. The result is shown in the skymap with a black 
cross as well as by a circle indicating its size. The best- 



fit coordinates are RAjaooo^SO^ 32"^ 20" ± ll|tat ±lllys 
and DECj2ooo=41° 30' 36'.'0 ± lf2stat ± I'^sys (for more 
details on the systematic uncertainties in the source po- 
sition determination, see Bretz et al. 2003). The posi- 
tion found is compatible within errors with the one de- 
termined by HEGRA, and barely compatible with the 
claims by Whipple mentioned above (in Konopelko et al. 
2007). 

The TeV J2032-t-4130 energy spectrum was obtained 
using the Tikhonov unfolding technique (Tikhonov & Ar- 
senin 1979). It can be fitted {x^/n.d.f = 0.3) by a power 
law function. The differential flux (TeV~^cm~^s~^) is: 

= (4.5 ± 0.3) X 10~^\E/1 Tey)-^ °±"-3 (1) 

The errors quoted are only statistical. The systematic er- 
ror is estimated to be 35% in the flux level and 0.2 in the 
photon index (see Albert et al. 2007a). The differential 
energy spectrum is shown in Figure [3l The HEGRA TeV 
J2032-H4130 and MAGIC Crab Nebula measured spectra 
(in Albert et al. 2007a) are shown in blue solid line and 
black dotted line, respectively. The MAGIC energy spec- 
trum is compatible both in flux level and photon index 
with the one measured by HEGRA. 

Crab Nebula data from the same periods and zenith 
angle distributions were studied with the same analy- 
sis chain to check for any systematic deviation due to 
the long observation period. No indication of time vari- 
ability was observed: the source integral flux is constant 
within errors, at 3% of the Crab Nebula flux. The rela- 
tive systematic uncertainty in the ratio of both fluxes was 
estimated to be less than 10%. This uncertainty comes 
mainly from the slightly different atmospheric transmis- 
sion conditions and differences in the detector parameters 
during data taking of the source and the Crab Nebula. 

For illustrative purposes, the dotted fines in Figure [3] 
represent one-zone hadronic and leptonic models of the 
high energy emission, both consistent with observations 
at lower energies in the region. Under the hadronic sce- 
nario, the tt'^ are obtained from a proton parent popula- 
tion described by a power law (index F = —2) with expo- 
nential cutoff at 100 TeV. The cutoff value was adopted 
to be consistent with the upper limit at the highest en- 
ergies coming from the HEGRA spectrum. The inverse 
Compton spectrum is obtained from an electron popula- 
tion with equal index and a 40 TeV exponential cutoff 
scattering off the CMB photons. As in Aharonian et 
al. 2005, we do not consider here the conditions under 
which particles are accelerated or how they lose energy. 
Our leptonic fits (see also the quoted paper for an SED 
representation) can only cope with the data if we are 
actually looking at a Compton peak around the energy 
range of detection, which is not fully discarded within er- 
rors. Both models are compatible with the high energy 
emission. Confirming the reality of the diffuse emission 
detected at lower energies is crucial to distinguish be- 
tween these and more complex models. 

4. CONCLUDING REMARKS 

MAGIC observations confirm the location of TeV 
J2032-h4130 found by HEGRA. The MAGIC observation 
shows an extended source with a significance of 5.6cr. We 
find a steady flux with no signiflcant variability within 
the three year span of the observations (with the flux be- 
ing at a similar level of the HEGRA data of the period 
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2002-2005). We also present the source energy spectrum 
obtained with the lowest energy threshold to date, which, 
within errors, is compatible with a single power law. 

We thank the lAC for the excellent working conditions 



at the ORM. The support of the German BMBF and 
MPG, the Itahan INFN, the Spanish CICYT, the ETH 
Research Grant TH 34/04 3, and the Polish MNil Grant 
1P03D01028 is gratefully acknowledged. 



REFERENCES 



Abdo, A.A. et al. 2007, ApJ, 658, 33. 

Aharonian F. A. et al. (HEGRA collaboration) 2002, A&A, 393, 
L37. 

Aharonian F. A. et al. (HEGRA collaboration) 2005, A&A, 431, 
197. 

Aharonian F. A. et al. (HESS collaboration) 2006, A&A, 460, 743. 

Albert J. et al. (MAGIC CoUab.), 2006a, Science 312, 1771. 

Albert J. et al. (MAGIC CoUab.), 2006b, ApJ, 637, L41. 

Albert J. et al. (MAGIC CoUab.), 2006c, ApJ, 643, L56. 

Albert J. et al. (MAGIC CoUab.), 2006d, ApJ, 638, LlOl. 

Albert, J. et al. (MAGIC Collaboration), 2007a, "VHE Gamma- 
Ray Observation of the Crab Nebula and Pulsar with MAGIC" 
accepted by ApJ, arXiv:0705.3244 

Albert J. et al. (MAGIC CoUab.) 2007b, " Very high energy 
gamma-ray observations during moonlight and twilight with 
the MAGIC telescope" submitted to Astroparticle Physics, 
astro-ph/0702475 

Albert J. et al. (MAGIC CoUab.) 2007c, "Implementation of 
the Random Forest Method for the Imaging Atmospheric 
Cherenkov Telescope MAGIC" submitted to Astroparticle 
Physics, arXiv:0709.3719 

Amenomori M. et al. 2006, Science, 314, 439. 

Anchordoqui, L. A. et al. 2007, 2007, PhRvD, 75, 063001. 

Bednarek, W. 2003, MNRAS, 345, 847. 

Bosch-Ramon V., Aharonian F. & Parodes J. M. 2005, A&A 432, 
609 



Bretz, T. & Wagner, W. (MAGIC Collaboration) 2003, Proc. of 

28th ICRC (Tsukuba), 2947. 
Butt Y. et al. 2003, ApJ 597, 494. 
Butt Y. et al. 2006, ApJ 643, 238. 

Cortina, J. et al. (MAGIC CoUab.), 2005, Proc. of the 29th ICRC, 

Pune, India, 5-359, astro-ph/0508274 
Domingo-Santaman'a, H. et al., 200S^ Froc. of 29th ICRC (Pune) 

5 363. 

Domingo-Santaman'a, E. & Torres, D. F. 2006, A&A 448, 613. 

Fegan, D. J., 1997, Journal of Physics G Nuclear Physics, 23, 1013. 

Hillas A. M., 1985, Proc. of the 19th ICRC, La JoUa, 3, 445. 

Horns D. et al. 2007, A&A 469, L17. 

Lang M. J. et al. 2004 A&A 423, 415. 

Marti J. et al. 2006, A&A 451, 1037. 

Marti J. et al. 2007, A&A 472, 557. 

Mukherjee R. et al. 2003, ApJ 589, 487. 

Konopelko A. et al. 2007, ApJ 658, 1062. 

Paredes J. M. et al. 2007 ApJ 654, L135. 

Rolke, W. A., Lopez, A. M.& Conrad, J., 2005, Nucl. Instrum. and 

Meth. A551, 493. 
Tikhonov A. N. & Arsenin V. Ja. Methods of Solution of Ill-posed 

Problem - M (Nauka, 1979). 
Torres D. F., Domingo-Santamarfa E., & Romero G. E. 2004, ApJ 

601, L75. 



MAGIC observations of the unidentified TeV 7-ray source TeV J2032+4130 



5 



TABLE 1 

Observing periods, zenith angle ranges and observation modes. 



Year 


T[h] 


Z.A.[dcg] 


Mode 


2005 


18.1 


13-30 


ON/OFF 


2006 


60.1 


11-44 


Wobble 


2007 


15.5 


11-30 


Wobble 
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TABLE 2 

Events rec;orded above SIZE> 800 pe. Won, N^ff and A^cx refer to the number of events recorded in the direction of the 

SOURCE, THE NORMALIZED BACKGROUND AND THE 7-RAY EXCESS RESPECTIVELY. THE NORMALIZATION RATIO /norm AND SIGNIFICANCE N„ 

ARE ALSO SHOWN. 



Yecir 




^off 


Nex 


/norm 




2005 


641 


576 


65 


0.47 


2.2 


2006 


688 


559 


129 


0.20 


4.8 


2007 


175 


136 


39 


0.20 


2.9 


Overall 


1504 


1271 


233 


0.27 


5.6 
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TABLE 3 

UL FOR SOURCES IN THE FOV, ABOVE 500 GeV. 



Target Name 


Flux (Crab Nebula) 


Cyg X-3 


0.011 


WR 146 


0.010 


3EG J2033+4118 


0.009 
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Fig. 1. — Distribution of the ^■^-parameter for events coming from the direction of TeV J2032+4130 (SIZE>800 pe), the background 
distribution subtracted (black points). A convolved radial Gaussian fit F=Axexp(-O.50^/((T^^y +cr^rc)) indicated by the solid black line 
with (Tsrc = 5.0±1.7 arcmin. The Tpsf was measured from MC simulation and validated with Crab Nebula observations to be iTps/=5.2±0.1 
arcmin (dashed black line). 
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Fig. 2. — Gaussian-smoothed ((t=4') map of 7-ray excess events (background-subtracted) for energies above 500 GeV. The MAGIC 
position is shown with a black cross. The surrounding black circle corresponds to the measured Icr width. The last position reported by 
Whipple is marked with a white cross while the HEGRA position is shown with a blue cross in the center of the field of view. The error 
bars, in all cases, correspond to the linear sum of the statistical and systematic errors. The green crosses correspond to the positions of 
Cyg X-3, WR 146 and the EGRET source 3EG J2033-I-4118. The ellipse around the EGRET source marks the 95% confidence contour. 
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Fig. 3. — Differential energy spectrum from TeV J2032+4130 as measured by the MAGIC telescope in black solid line. The grey shadow 
shows the Ic error in the fitted energy spectrum. The flux observed by Whipple in 2005 and in the Milagro scan are marked with colored 
squares (blue and pink, respectively). The grey dotted line represents the Crab Nebula energy spectrum measured by MAGIC. The blue 
line shows the HEGRA energy spectrum. Theoretical one-zone model predictions are depicted with dashed lines. 



